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We investigate the shear viscosity of leptons in a strongly coupled quark gluon plasma (QGP). 

We find that the shear viscosity due to the lepton-quark scattering is inversely proportional to the 
ratio of electric conductivity of the QGP to temperature up to the leading logarithmic order of the 
electromagnetic coupling. The finding implies that the thermal leptons form a more viscous fluid 
than the quarks by a large ratio. Using the known result for the electrical conductivity of strongly 
coupled plasmas obtained from gauge/gravity duality, we find that the lepton shear viscosity is sup¬ 
pressed compared with the one from lepton-lepton scattering. Gonsistently, we find an enhancement 
of the energy loss of hard leptons in a strongly coupled scenario compared with that in a weakly 
coupled plasma. 


The transport properties of thermal plasmas have been 
widely investigated for decades, because such plasmas 
play important roles in many areas of physics, rang¬ 
ing from thermonuclear fusion to astrophysics. Most re¬ 
cently, plasmas including free quarks have been studied in 
the context of relativistic heavy ion collisions and cosmol¬ 
ogy. In weakly coupled scenarios described by quantum 
electrodynamics (QED) or perturbative quantum chro¬ 
modynamics (QCD), transport coefficients such as the 
shear viscosity rj and electric conductivity a can be com¬ 
puted in the framework of kinetic theory Simi¬ 

lar studies have been also carried out for weakly cou¬ 
pled hadronic gas and a variety of quantum field theories 
[sl-fl^. In strongly coupled scenarios, the perturbative 
calculation and the quasi-particle description underlying 
kinetic theory become invalid. However, the transport 
coefficients can still be evaluated non-perturbatively us¬ 
ing Kubo formulas. For example, Demir and Bass m 
used the Kubo formalism combined with a phenomeno¬ 
logical transport model to compute the ratio of shear 
viscosity to entropy density (rj/s) of a thermal hadronic 
gas (see relevant studies of hadronic matter in [13, [II). 
The gauge/gravity duality can be employed to evalu¬ 
ate r] for strongly coupled plasmas via the Kubo for¬ 
mula BO giving rise to the well-known Kovtun-Son- 
Starinets (KSS) bound of ry/s. 

Here we address a more subtle situation - the mix¬ 
ture of both weakly coupled and strongly coupled sectors 
of the plasma - which exists in the practical cases such 
as the quark gluon plasma present in the early universe. 
The plasma comprises both leptons and colored quanta 
(quark and gluons), where the interactions among leptons 
are weakly coupled but the interactions among quarks 
and gluons are strongly coupled. The two sectors are 
connected by electroweak scattering between leptons and 
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quarks. Intuitively, one would expect a minor influence 
of the leptonic sector on the QCD sector, while the mag¬ 
nitude of the inverse effect is less obvious. We address 
this question by computing rj of light thermal leptons 
embedded in the strongly coupled plasma. 

Because of their weak coupling, the Boltzmann ap¬ 
proach can be applied to analyze the leptonic trans¬ 
port. There exist two collisional terms: one corre¬ 
sponds to lepton-lepton scattering and the other comes 
from lepton-quark scattering, which encodes the non- 
perturbative nature of the QCD sector. We denote the 
former contribution to the viscosity of the QED plasma 
by ?7qed and the latter by rjmix- Since the collisional 
terms are additive in the linear Boltzmann equation, one 
could compute rimix and rjQBD separately. The complete 
shear viscosity of leptons 

Vc ~ - (1) 

?7mix + T?QED 

will be always smaller than each individual contribution. 

We will now briefly describe our strategy to tackle the 
problem and mention our salient findings. In the Boltz¬ 
mann approach, the collisional term is proportional to 
the square of the scattering amplitude, which can be re¬ 
lated to the imaginary part of the photon self energy 
through the optical theorem. The photon self energy 
can be further written as the current-current correla¬ 
tion function, which can be computed via perturbative 
or non-perturbative approaches. The schematic figure 
is illustrated in FigUJ It was found in weak coupling 
that the collisional integral is dominated by an infrared 
(IR) divergence led by small momentum transfers. Such 
a divergence can be regularized in terms of a logarith¬ 
mic power of the coupling, ln(l/e) for QED and \n{l/g) 
for pQCD, by introducing an IR cutoff around the scale 
of Debye mass Qi, which gives the so-called leading 
logarithmic results of transport coefficients. 

In the strongly coupled scenario (in the QCD sector), 
we will use the correlation function obtained from the 
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FIG. 1: The self-energy diagram for thermal leptons. The 
curvy lines correspond to off-shell thermal photons and the 
solid lines represent the thermal leptons. The blob incorpo¬ 
rates the strong coupling in the QCD sector to all orders. 


gauge/gravity duality in the large number of colors Nc 
limit for the lepton-quark collision. Given that small- 
momentum transfers dominates we focus on the long- 
wavelength limit. It turns out that the collisional inte¬ 
gral is mainly contributed by collinear divergence when 
the intermediate virtual photons become almost lightlike. 
By ignoring the contribution from the longitudinal part 
of the correlator, we further relate rymix to the electric 
conductivity CTc of the color sector. We find that 

Vniix/T^ '^T/ac, ( 2 ) 


and the analogous equation for E<(p), where Sp and 
are Wightman functions of on-shell leptons and 
electromagnetic currents generated by thermal quarks 
in momentum space. In thermal equilibrium, one may 
rewrite in terms of the retarded correlator 

n^i/(9) linear response theory. Note that is 

proportional to since we only consider the leading- 
order interaction for the electromagnetic coupling. How¬ 
ever, contains all orders of the strong coupling g. 

One could easily show that r>(p) = express¬ 

ing detailed balance. 

The Boltzmann equation now becomes 


P^d^f{p,x) = 


tr 


e((7o)Im[n^^(g)] 

{27r)HEj,_g (g2)2 

i -b - ^) + m)Y] D{q,p) 


( 6 ) 


where Ep-q = + and 


Diq,p) = fip,x){l +f{q,x)){l-f{p-q,x)) 

f{p,x))f{q,x)f{p-q,x). (7) 


which is physically intuitive since the strength of colli¬ 
sions between two sectors is characterized by CTc in the 
long-wavelength limit. When a^T is large, rjmix will be 
suppressed making thermal leptons more viscous. 

We also study the energy loss for high-energy leptons 
or heavy leptons like muons, which slowly thermalize in 
the medium. We find that their scattering with the QCD 
sector results in a linear divergence in the IR regime for 
heavy leptons, which dominates over the logarithmic di¬ 
vergence from the scattering with thermal leptons. For 
light leptons, the 1/Nc suppression of ?7niix and the 0{Nc) 
enhancement of the energy loss mutually support each 
other. In the following, we now present the details of our 
approach and derivations. We denote the electric charge 
of the leptons by e and that of the single flavor of quarks 
by Q. The generalization to several quark flavors with 
differing electric charges is straightforward. 

In thermal equilibrium, the interaction rate of leptons 
in the relativistic Boltzmann approach can be written as 

^dpf{p,x) = -f{p,x)T>{p) -b (1 - f{p,x))T<{p), (3) 

H/p 


where /(p, x) denotes the distribution function of leptons 
in phase space. Here r>(p) and r<(p) represent the ra¬ 
diation and absorption rates, respectively. Based on the 
optical theorem, one finds 


r>«)(p) = (-)-_tr (|i-bTO)E>«)(p) 


(4) 


Here f{p,x) and f(p,x) correspond to fermionic and 
bosonic distributions and e(qo) denotes the sign function 
of go following the convention in [l^ . One can easily 
check that D{q,p) = 0 when f{p,x) and f{p,x) are the 
thermal distribution functions. 

Following the general approach in , we intro¬ 

duce the perturbations away from equilibrium for the lep¬ 
tons as / = h + 6f, where h(p, x) and n(p, x) correspond 
to the thermal distributions of fermions and bosons in 
the vicinity of the rest frame. We then assume that the 
perturbation takes the form, 

Sf{p, x) = (1 - n{p, x))n{p, x)xip, x), (8) 


where 


x(t,x) 


B{p) 

T 


p-fPdiUj 


(9) 


with Uj being the velocity and T being the temperature 
of the medium. Here B{p) is associated with g and the 
part corresponding to the bulk viscosity is dropped by 
setting V ■ u = 0 for simplicity. Collecting the leading- 
order contributions in (jS]), we find 



e(go) 


im [n^^(<?)] 


J {2Tr)3Ep.q[ 

n{q, x)h(p — q,x 
h{p, x) 


(g2)2 

q) 




,( 10 ) 


where m is the lepton mass. The are given by 

T,>(p) = j q)j'', (5) 


where 

Bij{p,q) 


(2p'^(p - g)" - p ■ (p - + mV"), 

{B,j{p) - Bij{p - q)) , 
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B,j{p) = B{p) . (11) 

From the definition of 77 as the perturbation of the energy- 
stress tensor, one obtains 


correlators [3^. Now, we may choose the simple ansatz 

lr,|2 


Bip) = C 


r2 In(^) 


(17) 


The primary task now is to solve the integral equation 
for B{p). In general, the integral equation must be 
solved numerically m- Nonetheless, in order to make 
manifest comparison with the result of the QED plasma, 
we may choose an appropriate ansatz for B{p) and work 
out an approximate analytic solution of Tymix up to the 
leading logarithmic order. 

From the structure of the integrand in the momentum 
integral of (nni), one may anticipate the dominance of low 
momentum transfer. This is true for QED and weakly 
coupled QCD plasmas. It also turns out to be the case 
when the QCD sector is strongly coupled. Therefore, we 
will work in the long-wavelength region for light leptons 
(m < T), 


go < |q| < T < |p|. (13) 

Eor convenience, we may further decompose into 

the transverse and longitudinal parts, 

(14) 

where are projection operators. In the long- 

wavelength limit, (nni) reduces to 


( 


3 



d|q|(is 

27r2|q|3s 


' r|p|(l-s^)e(go) 


Sy(p,q)Im [n^(g) - n^(g)] 


,(15) 


where we preserve the small lepton mass in to avoid 
the collinear divergence at s = cos 0 = ±1 in the photon 
propagators [s^- 

In order to study the effect of the non-perturbative dy¬ 
namics in the QCD sector, we take the D3/D7-brane sys - 
tem in gauge/gravity duality as a concrete example | 1 § |. 
which is dual to a strongly coupled supersymmetric gauge 
plasma with quarks in the fundamental representation. 
Such a system as an analogue of the strongly coupled 
QGP in holography was frequently used to study elec- 
tro mag netic signatures emitted from a thermal medium 
[l9l - l^ . In the long-wavelength limit for q = (go, 0,0, g^), 
one finds 


where C = 2^^/{e^QpNf) is a dimensionless constant, 
which serves as the leading-order solution of (HHD in the 
long-wavelength limit. 

B (p) diverges in the limit m —^ 0 owing to the collinear 
divergence in the integral over s. As a regularization we 
may take m ^ eT as the thermal mass of leptons, which 
is equivalent to introducing an IR cutoff at |q| ^ eT to 
regularize the IR divergence in QED. We could further 
take ln( 2 |p|/(eT)) « ln(l/e) as an approximation ap¬ 
plied in QED up to the leading logarithmic order. There 
exists a caveat that such an approximation may break 
down when |p| > T/e. For quantitative correctness, one 
may also have to incorporate the Landau-Pomeranchuk- 
Migdal (LPM) effect with multiple Coulomb scattering 
for high-energy leptons. 

Inserting (II3 into we obtain the Tymix in a 

strongly coupled supersymmetric plasma modeled by the 
D3/D7-brane system, 

D 3 /D 7 ^ _ yio'i 

Vix ^ 30iV^Ar^e2Q2in(I/e)’ ^ > 

where Ip ~ 116. To make a comparison with tjqed on 
equal footing, one could follow the same approach and 
utilize the photon self energy in O(e^) to solve the Boltz¬ 
mann equation with a similar ansatz, which gives rise to 

M 


= 5nenn{l/ey 

Comparing m with (HID we conclude that the ratio 
^mfxis suppressed in the large-A"c limit. 

Based on the dominance of the collinear divergence 
in our derivation, we can generalize our result by con¬ 
necting Im[n^(g)] in the long-wavelength limit with the 
direct-current (DC) electrical conductivity Uc of the QCD 
sector. The collinear divergence occurs when go ^ |q|, 
which yields Im[n^(g)] —>■ 0. Therefore, the integral on 
the left hand side of (fT5l) will be dominated by the contri¬ 
bution from Im[n^(g)]. Recall that the DC conductivity 
could be dehned as 


o-c 


1 £(go)x(go,|q|)|qhgo 

4 go->0 go 


( 20 ) 


where x(go, |q|) denotes the spectral density, which can 
be computed from the retarded correlator. 


Im [n^(g)-n^(g)] 


Q^N^NfT g 2 e(go) 
47r go 


where the extra factor e(go) is introduced for the cut 


x{qo, |q|) = - 2 im ( 2 n^(g) -h n^(g)). (21) 

The static limit (go —>■ 0) is approximately equivalent to 
the long-wavelength limit here. By dropping Im[n^(g)] 














4 


and replacing Im[n^(( 7 )] with the DC conductivity as 


^ / e(go)Im[n^(g)] \ 

we thus obtain the more general result: 

T^Ip 

??mix l20(Tce^ ln(l/e) 


( 22 ) 


(23) 


’ll 

Hqed 



Armed with (1^51) . we can evaluate in the lead¬ 
ing logarithmic order in different systems provided CTc is 
known. In pQCD at finite temperature, the DC conduc¬ 
tivity up to the leading logarithmic order mil leads 
to 


pQCD ^ T^Ip{N^-l)gHn{l/g) 

^ coN^Nfe^QHn{l/e) ’ ^ ’ 

where Cq is a numerical constant [s^. By contrast, in the 
Sakai-Sugimoto (SS) model of holographic QCD 
the known value of the DC conductivity in the deconfined 
phase [28l - [30l| yields 


'/mix 


QttIpT'^Mkk 
AOXN^ e^Q^ ln(l/e)’ 


(25) 


where A denotes the t’Hooft coupling and Mkk corre¬ 
sponds to the Kaluza-Klein mass as the meson mass scale 
in the gauge theory. Qualitatively, the ratio rymix/ ?7qed in 
pQCD and that in the SS model despite the competition 
of Mkk and T are g'^hi{\/g) and l/(AtA^c) suppressed, 
respectively. These results may be representative of the 
viscous behavior of thermal leptons in high-temperature 
(pQCD) and strongly coupled, intermediate-temperature 
regimes of QGP. 

Using the analytical expressions derived above, we 
may attempt a crude estimate of the ratio rymix in the 
cosmological QGP beyond the deconfinement transition 
temperature. According to various approaches for the 
electrical conductivity of QGP [3l|, one finds CTc/T ^ 
0.02 - 0.1 from T = 0.2 GeV to T = 0.6 GeV. This 
implies a ratio 77mix/?7QED ^ 0.12 — 0.6, corresponding 
to rjc/VQED ^ 0.11 — 0.38. We conclude that the lep- 
tonic visosity in cosmological QGP was dominated by the 
interactions between leptons and quarks. Nevertheless, 
UcjT drops rapidly near the deconfinement transition, 
where the lepton-lepton scattering suppresses the lepton- 
quark interaction and rjcIVQEB thus increases. As shown 
in Fig|2 by utilizing OcjT obtained from the stro ngly 
coupled non-conformal gauge theory in holography [32| . 
Vc/VQEi) rapidly increases near the transition tempera¬ 
ture. 

We now apply the same considerations to the interac¬ 
tion of hard leptons with the QGP. Because here we do 
not assume that the leptons are thermal, this calculation 
also has relevance to the QGP produced in relativistic 
heavy ion collisions. For energetic or heavy leptons such 
as muons which slowly thermalize, we may compute the 


FIG. 2: The ratio of the complete shear viscosity of leptons in 
a strongly coupled non-conformal plasma to that in the QED 
plasma versus temperature. 


radiation rate T^{Ep) by simply replacing the thermal 
cut- propagator of leptons in the previous computations 
with the one in vacuum. Following the definition in [l^ 
with the long-wavelength approximation, the energy loss 
for a hard lepton traveling along the z direction in the 
strongly coupled supersymmetric plasma characterized 
by the D3/D7 system could be written as 


dEp 

dz 


e^Q^NcNfT 

IbTT^ 


d|q|c?s 


(l -I- D^lqp) 

+ D^lqp 


/ u^(l — s^) v{l — v^) \ 

l^(l_„2g2)2 + (l_„2g2) ) ’ 


(26) 


where s = cosd and v = \p\/Ep and D = (27rT)“^ is the 
diffusion constant. We focus on the large-velocity and 
small-velocity regions, where the analytic expressions are 
accessible. For light leptons m ^ eT <C |p|, one could 
simply drop the diffusion terms and obtain 


/dEp\ 

\ /D3/D7 


\n{l/e)N^NfT^ 
Stt^ 


for V —>■ 1(27) 


For heavy leptons with v <C D|q|, the presence of the 
diffusion term prevents the 1 /v divergence in (1261) , while 
there still exists linear divergence in the IR regime for the 
|q| integral. By simply taking the IR cutoff qm eT, 
we acquire 


fdEp\ eQ^N.NfT^ ^ T 

~-;;—ft— for r' —. 

\ dz J dirEp Ep 


Compared with the energy loss in the leading order of 
two limits in the QED plasma. 


dEp^ 

) ^ 

e'^ln(l/e)T^ 

for 

u —>■ 1, 

dz y 

' QED 

2477 



dEp'' 

) ^ 

e'‘ln(l/e)T3 

for 

T 

V —>• — 

dz y 

' QED 

SGnEp 


Ep 


it turns out that the pattern of the energy loss of light 
leptons in a QED plasma is similar to that in the strongly 
coupled supersymmetric plasma. Nonetheless, despite 
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the prefactor altered by the choice of qm, the diffusion 
at small momentum transfer further enhances the energy 
loss for heavy leptons in the strongly coupled scenario. 
Even for light leptons, the energy loss is 0{Nc) enhanced 
in the D3/D7 system, which is in accordance with the 
l/Nc suppression of 

In summary, we have studied the shear viscosity and 
energy loss of leptons in strongly coupled (weakly cou¬ 
pled) QCD-like plasmas. It turns out that the non- 
perturbative effect of the QCD sector substantially af¬ 
fects the leptonic transport, rendering the lepton fluid 
more viscous. In order to obtain more accurate results 
for the leptonic viscosity, numerical solvers of the Boltz¬ 
mann equation beyond the leading logarithmic order are 
required. In addition, more realistic simulations of the 
current-current correlator should be involved. Our ap¬ 
proach can be easily generalized to other fluids that 
weakly interact with a strongly coupled medium. For in¬ 
stance, in the semi-holographic model of QGP [sj, where 
hard and soft gluons are connected by effective pertur¬ 
bative couplings, our approach could be utilized to study 


the transport of hard gluons within the QGP. Moreover, 
there exist also thermal photons in the cosmic plasma, 
the study of the non-perturbative effect from the QGD 
sector on thermal-photon transport in the similar frame¬ 
work could be an interesting issue. 
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